		twoway	(scatter $p vars1 if dataset == "bs_decomposition_SS_ProdQ_new_klme_TFPQIV", mcolor($cl) msymbol($sym)) ///
				(rcap lows_$p highs_$p vars1 if id == 2 & dataset == "bs_decomposition_SS_ProdQ_new_klme_TFPQIV", lwidth(medthick) lcolor($cl ) color(%50)) ///
				(rcap lows_$p highs_$p vars1 if id == 1 & dataset == "bs_decomposition_SS_ProdQ_new_klme_TFPQIV", lwidth(thin) color(gs6) ) ///
				(scatter $p vars2 if dataset == "bs_decomposition_SS_ProdQ_new_klme_TFPQIV_2lags", mcolor($cl) msymbol(Oh)) ///
				(rcap lows_$p highs_$p vars2 if id == 2 & dataset == "bs_decomposition_SS_ProdQ_new_klme_TFPQIV_2lags", lwidth(medthick) lcolor($cl ) color(%50)) ///
				(rcap lows_$p highs_$p vars2 if id == 1 & dataset == "bs_decomposition_SS_ProdQ_new_klme_TFPQIV_2lags", lwidth(thin) color(gs6) ) ///
				(scatter $p vars3 if dataset == "bs_decomposition_SS_ProdQ_new_klme_TFPQIV_rshare", mcolor($cl) msymbol(Dh)) ///
				(rcap lows_$p highs_$p vars3 if id == 2 & dataset == "bs_decomposition_SS_ProdQ_new_klme_TFPQIV_rshare", lwidth(medthick) lcolor($cl ) color(%50)) ///
				(rcap lows_$p highs_$p vars3 if id == 1 & dataset == "bs_decomposition_SS_ProdQ_new_klme_TFPQIV_rshare", lwidth(thin) color(gs6) ) ///
				(scatter $p vars4 if dataset == "bs_decomposition_SS_CS_klme_nic4", mcolor($cl) msymbol(Th)) ///
				(rcap lows_$p highs_$p vars4 if id == 2 & dataset == "bs_decomposition_SS_CS_klme_nic4", lwidth(medthick) lcolor($cl ) color(%50)) ///
				(rcap lows_$p highs_$p vars4 if id == 1 & dataset == "bs_decomposition_SS_CS_klme_nic4", lwidth(thin) color(gs6) ) ///
				(scatter $p vars5 if dataset == "bs_decomposition_SS_CS_klme_nic2_by_year", mcolor($cl) msymbol(Sh)) ///
				(rcap lows_$p highs_$p vars5 if id == 2 & dataset == "bs_decomposition_SS_CS_klme_nic2_by_year", lwidth(medthick) lcolor($cl ) color(%50)) ///
				(rcap lows_$p highs_$p vars5 if id == 1 & dataset == "bs_decomposition_SS_CS_klme_nic2_by_year", lwidth(thin) color(gs6) ) ///
				(scatter $p vars6 if dataset == "bs_decomposition_SS_ProdQ_new_klm", mcolor($cl) msymbol(Sh)) ///
				(rcap lows_$p highs_$p vars6 if id == 2 & dataset == "bs_decomposition_SS_ProdQ_new_klm", lwidth(medthick) lcolor($cl ) color(%50)) ///
				(rcap lows_$p highs_$p vars6 if id == 1 & dataset == "bs_decomposition_SS_ProdQ_new_klm", lwidth(thin) color(gs6) ) ///
				, xlabel(2 "R. TFP" 4 "M. Demand" 6 "M. TFPQ" 8 "Covariance" 10 "Scope" 12 "Adding" 14 "Dropping", angle(35)) ///
				ylabel(,grid) ///
				legend(order(1 "KLME" 4 "KLME 2Lags" 7 "KLME RSHARE" 10 "CS KLME NIC4" 13 "CS KLME NIC2 by Year" 16 "Q KLM") size(medlarge) rows(3) region(lcolor(white)) position(6) ring(1)) scheme(s1color) 
				*title("GO $p") 

			graph export "$figures/robustness_GO_$p.png", as(png) replace width(1800) height(1000)
